Optical absorption from a non degenerate polaron gas 
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We apply the Dynamical Mean Field Theory to calculate the low-density limit of the optical 
conductivity in the Holstein model. This non perturbative treatment allows to span continuously 
from the weak-coupling quasi-free electron behaviour to the strong-coupling small polaron picture. 
New features arise in the intermediate coupling regime: i) a noticeable transfer of spectral weight 
towards low frequencies occurs at low temperatures, when coherent processes acquire importance; ii) 
thermally activated Polaron Interband Transitions appear at non zero temperature, whose spectral 
weight has a non monotonic temperature behaviour. We propose an interpretation of the infrared 
active vibrations (IRAV) observed in the optical spectra of the underdoped superconducting cuprates 
in terms of small polarons in the intermediate coupling regime. 



The calculation of the optical conductivity of small po- 
larons is an important and unsolved problem in solid 
state physics. Although many approaches exist which 
are appropriate in specific limiting cases, an analytical 
treatment which is able to give the correct physical re- 
sults in all the regimes of parameters is still lacking. The 
simplest model which describes small polarons is the Hol- 
stein Molecular Crystal Model Q| , where the carriers in- 
teract with phonon modes localized on each lattice cell. 
The Holstein model was recently solved [|[ in the single- 
electron case by Dynamical Mean Field Theory (DMFT) 
[^ , which becomes exact in the limit of infinite lattice co- 
ordination. When applied to real lattices, this approach 
can be viewed as a valid interpolation scheme, provided 
that an appropriate finite dimensional density of states 
(DOS) is chosen ^^. The solution reported in ref. ||] 
is essentially analytical, and it overcomes the well-known 
difficulties inherent to numerical methods. 

The aim of this work is to generalize the exact DMFT 
solution of the Holstein model to the calculation of the 
optical conductivity. We shall focus on the most inter- 
esting intermediate coupling regime near the adiabatic 
limit, where novel features have been found in the sin- 
gle particle propagator, so that new interesting proper- 
ties can be expected in the optical spectra. Contrary 
to the usual strong coupling situation |p| , where the po- 
laron absorption is concentrated in a "band" located at 
frequencies much higher than the phonon frequency, we 
will show that at intermediate coupling strengths, narrow 
absorption peaks appear close to the phonon frequencies, 
which are due to transitions between different subbands 
in the polaron excitation spectrum. The presence of such 
low-frequency features is a very general property of finite 
band electron-phonon (e-ph) systems in the intermediate 



coupling regime, which could help interpreting the far 
infra-red spectra of the superconducting cuprates at low 
doping levels, as discussed in the conclusion. 

In the Holstein hamiltonian |l| , tight-binding electrons 
with hopping amplitude t are coupled locally to Einstein 
phonons with energy ujq • 
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At zero temperature, a crossover occurs from a quasi free 
electron to a small polaron as the e-ph coupling constant 
g increases. This crossover depends strongly on the adia- 
baticity parameter ujo/t [Q||]. If we restrict to the nearly 
adiabatic regime {ojo/t < 1), the crossover takes place 
when the polaron energy Ep = g'^ /ujq is of the order of 
the free electron bandwidth, i.e. for A = g'^/ujQt ~ 1. 
Although these arguments are based on the ground-state 
properties alone, it was pointed out in reference that 
the mechanism of polaron formation can be better un- 
derstood by analyzing the whole excitation spectrum, as 
deduced from the knowledge of the interacting propaga- 
tor G. It was shown that the polaron crossover is char- 
acterized by the successive opening of energy gaps in the 
low energy part of the spectral density as the e-ph cou- 
pling is increased. At intermediate values of A, although 
the polaron mass is still close to the free electron value, 
the spectral density starts to exhibit typical polaronic 
features, i.e. narrow subbands at low energies, together 
with a broad contribution at higher energies. The width 
W of the first subband, which gives a measure of the 
coherent (band-like) behaviour of the polaron, is rapidly 



reduced as A moves towards the strong coupling limit. 
This picture of the polaron crossover has important con- 
sequences on the optical properties. 

The calculation of the optical conductivity within the 
DMFT is simplified since, for symmetry reasons and due 
to the locality of vertices, no current vertex corrections 
enter in the Kubo formula ||9|,[l0|. Writing the electron 
spectral function as 
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then the optical conductivity per particle takes the form 

Cf{uj) = deN^(j)e Idv {fv- fuj+u)Pe{v)Pc{v + LU) (3) 



(2/7r)\/t2 - e2 is the free DOS of the real 
{t^ — e2)/3 is the corresponding current ver- 



where N^^ 
lattice, (/>£ 

tex ||ll| (which ensures a sum rule for the total spectral 
weight [|2|), / is the Fermi function, n is the particle 
density and ^ is a numerical constant that we shall drop 
in the following discussion M]. Taking the low density 
limit [n — > 0) at finite temperatures, following the pro- 
cedure described in ref. [H, eq. (0) becomes 



aH = -(1 



where we have defined 



^ AA(/3) 
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AA(/3)= j deN, f due-^^"-^"^ pi°\i^) (6) 



and p^^' is the spectral function at zero density (these 
quantities have been rescaled by e'^^" to obtain a finite 
result in the low temperature limit). The integrals in e 
can be performed analytically, while the remaining inte- 
grals in v are computed numerically, once the self-energy 
S is obtained as a continued fraction expansion by ap- 
plying the iteration scheme reported in reference g] . 

Our treatment applies to a dilute non- degenerate po- 
laron gas, where exchange effects can be neglected. Turn- 
ing to real systems, we expect that our results will be rel- 
evant at low enough density, i.e. when the Fermi temper- 
ature is negligible compared to the important tempera- 
ture scales of the problem (e.g. the temperature hujo/kB, 
below which phonon quantum effects become important 

Figure Q^ shows the results obtained at cja/t = 0.1 and 
A = 0.9, i.e. close to the polaron crossover |13]. The opti- 
cal response shows a typical pattern of subbands roughly 
spaced by ujq, which tends to become smooth at higher 
frequencies. This behavior is a characteristic signature 
of the polaron crossover M . At temperatures well above 



ujQ, the results are in good agreement with Reik's for- 
mula PI, showing that the physics in this limit is domi- 
nated by activated hopping processes. As the tempera- 
ture is decreased, a substantial transfer of spectral weight 
takes place, the low-frequency features becoming more 
pronounced. 




FIG. 1. Optical conductivity per particle at A = 0.9 and 
uo/t — 0.1, at different temperatures. Reik's result is shown 
at T = 5ujQ for comparison. Note the transfer of spectral 
weight occurring at low T (for the definition of Tw, see text). 
The curve marked T « Tw has been obtained with eq. (M) 

To get some insight on this peculiar behavior, it is 
convenient to analyze the polaron excitation spectrum 
by separating explicitly the contributions p^ and p^ at 
energies respectively above and below _Eo. At low tem- 
peratures, p^ has quasiparticle poles at energies below 
Eq + ujq plus a more or less incoherent background due 
to phonon scattering at higher energies. On the other 
hand, p^ represents incoherent states where one or more 
phonons are taken from the thermal bath. Correspond- 
ingly, by looking at the definitions (||) and (^), one can 
define A/" = A/'>+A/'< imdV = V» +V<> +V« . The 
term !?>> has contributions associated to both coherent 
band motion and phonon emission processes, while T)^^ 
and !?<< account for the thermally activated hopping. 

An estimate for the temperature Ty/ below which co- 
herent processes acquire importance can be obtained by 
analyzing the relative weight of N^ and M^ . Due to 
the exponential factor in eq. (pi), we see that at low 
temperatures N^ weights the coherent states close to 
Eq. Using the semicircular free DOS, we obtain N^ = 
yj2(m* jm)l-K (r/t)'^/^, where m* /m is the ratio of the 
polaron mass to the free electron value. On the other 
hand, being associated to incoherent states, the spectral 
function pf{v) is essentially e- independent and peaked 
around a discrete set of energies v = Eq + W — ploq, with 
a weight proportional to the thermal occupation factor 
g-P/3ajo ^ (recall that W is the width of the first sub- 
band in the spectral density, and p = 1,2,.. is the num- 
ber of thermal phonons). When evaluating the integral 
in A/"^ , such vanishing occupation factors are partly com- 



which gives rise 



pensated by the exponential e~^^'^~^°\ 
to an overall behavior A/"^ oc e"'"^. 

We find that the crossover between the two regimes 
occurs at T\y ~ W/log{t/W). This is not surprising, 
since as was noted in the introduction, W is the typical 
energy scale for coherent excitations (in the particular 
case studied here, where m* /m ~ 1.5 and W = O.To^o, 
we get Tw ^ O.Iwq). For T < Tw, only the terms V» 
and A/"^ survive in formula (0), which becomes: 
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where the spectral function p is evaluated at e = e{q = 0). 
In this regime, the conductivity scales as T/uj for any 
w > T, so that most of the spectral weight is transferred 
to the Drude peak a{uj) ~ S{uj). 

At higher temperatures, the main contribution to the 
optical absorption involves phonons that are already 
present in the thermal bath, so that one expects some 
spectral weight to move towards higher frequencies as the 
temperature is raised above T-^/. This is indeed the case, 
as can be seen in figure 1. Let us stress that this crossover 
from coherent to activated conductivity is a signature of 
the intermediate coupling regime. In the strong coupling 
limit, where W —> 0, the crossover temperature vanishes 
and phonon-assisted contributions always dominate. 
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FIG. 2. The same as preceding figure but in the region close 
to ojo and with a finer temperature scan. The curve T <^ Tw 
corresponds to T = 0.025a;o (the Drude peak is not shown). 

In Fig. we have reported the optical conductivity in 
the frequency range around loq, at temperatures close to 
Tw- Despite its apparent complexity, the observed pat- 
tern can be qualitatively understood by looking at the 
excitation spectrum close to the ground state, which is 
sketched in figure 3 versus the free-particle energy e (see 
eq. (0)). The dark curves are the quasiparticle states in 
the lowest order polaron subbands close to the ground- 
state, while the region above uq represents the incoher- 
ent background due to scattering from Einstein phonons 
(the states below Eq, which give rise to the hopping- 
like behavior at T ^ Tw, are not shown in fig. 3). The 



transitions from the ground-state to the incoherent back- 
ground, denoted by (a), give rise to a continuum of ab- 
sorption above wq- A transition from the ground-state to 
the bottom of the second subband is also possible (b), but 
this is hardly visible in fig. 2 owing to its negligible spec- 
tral weight. In addition, the Drude peak corresponding 
to coherent band motion is clearly visible in fig. 2. 
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FIG. 3. The dispersion of the lowest order polaronic sub- 
bands versus the free particle energy e (parameters are the 
same as in the preceding figures, and T = 0.2(i;o). The points 
are distributed according to the spectral function p^iy). 



The most striking feature in figure 2 is probably the 
appearance of a resonance below luq (c) , which is due to 
Polaron Interhand Transitions (PIT), as can be seen in 
figure 3. Indeed, because of thermal activation, long-lived 
states at energy Eq < v < Eq + ujq can be occupied with 
a probability e~^^'^~^^\ According to the Fermi golden 
rule, the transition probability at a given frequency uj is 
proportional to the joint density of the initial and final 
states separated by u. If we focus on transitions be- 
tween different subbands, this function has a Van Hove 
singularity at that value oj* where the band dispersions 
are parallel, which identifies an interband threshold (see 
arrow (c)). This resonance has a non monotonic temper- 
ature dependence, which can be understood as follows: 
the occupation probability of the initial state gives rise 
to an activated behavior at low T; as the temperature is 
raised above uj* , the difference in population of the initial 
and final states is balanced, which depresses the optical 
absorption according to the prefactor '--^ (1 — e"^'^ ) in 
formula (0). According to the discussion above, we ar- 
gue that PIT are not peculiar to the Holstein model, 
nor are they a consequence of the approximations used. 
Rather, they are closely related to the presence of several 
subbands in the spectral density at low energy, which 
should be a common feature of any finite band model 
in the intermediate e-ph coupling regime. Such inter- 
band transitions are obviously absent at weak coupling, 
where there are no subbands in the excitation spectrum, 
while in strong coupling they can hardly be resolved from 
the incoherent multiphonon absorption, due to their rel- 



atively small spectral weight. On the other hand, the 
results are not qualitatively affected by variations of the 
adiabaticity parameter, provided that uja/t < 1 ||]. 

In this work, we have analyzed the optical conductivity 
of a dilute gas of small polarons in the framework of the 
DMFT, focusing on the crossover regime between quasi- 
free electrons and small polarons which takes place at 
intermediate values of the e-ph coupling strength. The 
spectra at high temperatures do not differ significantly 
from the classical picture of hopping-likc conductivity, 
which predicts a featureless absorption band at finite fre- 
quencies. However, new interesting results appear at low 
temperatures, where the physics is governed by two dif- 
ferent energy scales. At temperatures below ujq, phonon 
quantum effects become important, and the optical spec- 
tra start to exhibit narrow peaks at low frequencies, to- 
gether with a broad background at higher frequencies. 
Upon further decreasing the temperature, a noticeable 
transfer of spectral weight takes place towards low fre- 
quencies. This occurs at a temperature Tw which is re- 
lated to W , the width of the first polaron subband in the 
spectral density, i.e. the typical energy scale governing 
coherent behaviour. 

Our results in the low frequency region (see Fig. 2) 
bear a strong resemblance with the far infrared spectra 
of both hole- and electron-doped cuprates at low dop- 
ing levels, where narrow peaks are observed at typical 
phonon frequencies, usually called infrared active vibra- 
tions (IRAV) |0,0]. From a theoretical point of view, 
if we focus on such low energy features, i.e. much lower 
than the typical energies of electronic correlations (e.g. 
J in the t-J model), the inclusion of correlation effects 
would only enter through a renormalization of the hop- 
ping parameter. Therefore, we expect that a description 
in terms of the Holstein model with renormalized pa- 
rameters should correctly account for all the observed 
features in the frequency range of interest |19[| . On the 
other hand, the present single-polaron treatment, which 
neglects polaron-polaron interactions, is appropriate at 
low doping levels, i.e. precisely where the IRAV are best 
resolved. 

In the measured optical spectra, the IRAV peaks ap- 
pear as shoulders on both sides of the main phonon lines, 
and they have been ascribed to the lowest order tran- 
sitions inside the polaron potential- well [g^l- We have 
shown here that in the intermediate e-ph coupling regime, 
the polaron absorption exhibits narrow peaks located 
above and below the phonon frequency |2l| , whose T de- 
pendence is non monotonic. These peaks come from two 
different processes: transitions to incoherent scattering 
states (denoted by (a)), and thermally activated Polaron 
Interband Transitions (denoted by (c)). Although indica- 
tions of a non monotonic behavior of some IRAV already 
exist |18| , p2[ , a more detailed experimental analysis of 
their temperature dependence is needed to confirm these 
ideas. 
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